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The  coherence  function  of  sound  waves  propagating  through  an  intermittently  turbulent  atmosphere 
is  calculated  theoretically.  Intermittency  mechanisms  due  to  both  the  turbulent  energy  cascade 
(intrinsic  intermittency)  and  spatially  uneven  production  (global  intermittency)  are  modeled  using 
ensembles  of  quasiwavelets  (QWs),  which  are  analogous  to  turbulent  eddies.  The  intrinsic 
intermittency  is  associated  with  decreasing  spatial  density  (packing  fraction)  of  the  QWs  with 
decreasing  size.  Global  intermittency  is  introduced  by  allowing  the  local  strength  of  the  turbulence, 
as  manifested  by  the  amplitudes  of  the  QWs,  to  vary  in  space  according  to  superimposed  Markov 
processes.  The  resulting  turbulence  spectrum  is  then  used  to  evaluate  the  coherence  function  of  a 
plane  sound  wave  undergoing  line-of-sight  propagation.  Predictions  are  made  by  a  general 
simulation  method  and  by  an  analytical  derivation  valid  in  the  limit  of  Gaussian  fluctuations  in 
signal  phase.  It  is  shown  that  the  average  coherence  function  increases  as  a  result  of  both  intrinsic 
and  global  intermittency.  When  global  intermittency  is  very  strong,  signal  phase  fluctuations  become 
highly  non-Gaussian  and  the  average  coherence  is  dominated  by  episodes  with  weak  turbulence. 
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I.  INTRODUCTION 

The  main  purpose  of  this  article  is  to  examine  how  tur¬ 
bulent  intermittency  affects  the  coherence  of  propagating 
sound  waves.  Coherence  of  sound  in  a  turbulent  atmosphere 
is  important  for  remote  sensing  1  and  assessment  of  perfor¬ 
mance  of  acoustic  arrays  for  source  localization."  Theoretical 
treatments  of  the  coherence  function  (which  describes  the 
dependence  of  coherence  on  spatial  separation  between  the 
observation  points)  for  waves  propagating  in  a  random  me¬ 
dium  are  well  developed  in  many  regards. 3-6  The  calculated 
coherence  function  depends  significantly  on  the  model  for 
the  turbulence.  In  recent  years,  turbulence  models  have  be¬ 
come  more  realistic  in  contrast  to  early  studies  that  dealt 
with  homogeneous,  isotropic  random  scalar  fields  possessing 
a  Gaussian  correlation  function.  For  example,  the  coherence 
function  has  now  been  calculated  for  the  case  of  inhomoge¬ 
neous,  anisotropic  turbulence  with  realistic  spectra  of  tem- 
perature  and  wind  velocity  fluctuations.  ’  Here,  we  consider 
the  introduction  of  an  additional  realism,  namely,  intermit¬ 
tency. 

'Portions  of  this  work  were  presented  in  V.  E.  Ostashev,  D.  K.  Wilson,  and 
G.  H.  Goedecke  “Intermittent  scalar  QW  model  and  sound  propagation 
through  intermittent  turbulence,"  in  Proceedings  of  the  12th  Long  Range 
Sound  Propagation  Symposium,  New  Orleans,  LA,  2006,  pp.  429-442. 
^Electronic  mail:  d.keith.wilson@usace.army.mil 
Electronic  mail:  vladimir.ostashev@noaa.gov 
d  Electronic  mail:  ggoedeck@nmsu.edu 


Intermittency  refers  to  the  tendency  of  turbulence  to  oc¬ 
cur  in  spatial  and  temporal  bursts  of  activity.  This  property  of 
turbulence9'10  plays  an  important  role  in  many  practical  prob¬ 
lems.  Mahrt11  has  proposed  classifying  intermittency  as  in¬ 
trinsic  or  global.  The  former  occurs  on  scales  less  than  the 
outer  scale  of  turbulence  (the  scale  of  largest  eddies)  as  the 
turbulent  cascade  process  progressively  concentrates  turbu- 
lent  energy  dissipation  into  smaller  regions  of  space.  ‘  Glo¬ 
bal  intermittency  occurs  on  scales  larger  than  the  outer  scale 
and  may  have  several  causes  related  to  uneven  production  of 
turbulence  in  a  particular  environment.  Possible  causes  in¬ 
clude  wind  gusts  from  large-scale  convective  systems  or  to¬ 
pographic  flow,  irregular  episodes  of  turbulent  mixing  in  sta¬ 
bly  stratified  (night  time)  boundary  layers,  large  organized 
coherent  structures,  and  uneven  heating  of  the  ground  due  to 
clouds  and  variations  in  ground-surface  properties.  The  terms 
small  scale  and  large  scale  are  also  used  to  describe  intrinsic 

11  13 

and  global  intermittency,  respectively.  ' 

Figure  1  is  a  conceptual  illustration  of  intermittency  and 
its  effects  on  wave  propagation.  The  eddies  occur  in 
“clouds,”  each  of  which  represents  a  global  intermittency 
event.  Each  turbulent  cloud  can  be  regarded  as  having  its 
own  outer  scale  and  turbulent  strength.  Within  each  cloud, 
the  smaller  eddies  organize  into  intermittent  patches  as  a 
result  of  the  turbulent  cascade  process  (intrinsic  intermit¬ 
tency).  The  scattering  experienced  by  sound  waves  propagat¬ 
ing  along  various  paths  through  this  region  varies  greatly. 


J.  Acoust.  Soc.  Am.  124  (2),  August  2008 


0001  -4966/2008/1 24(2)/743/1 5/$23.00 


©  2008  Acoustical  Society  of  America  743 


FIG.  1.  (Color  online)  Illustration  of  intrinsic  and  global  intermittency,  and 
their  effects  on  wave  propagation.  The  solid  circles  represent  turbulent  ed¬ 
dies  and  the  dashed  lines  represent  global  regions  of  active  turbulence.  The 
numbers  and  arrows  indicate  illustrative  propagation  paths  as  described  in 
the  text. 


according  to  the  number  and  intensity  of  clouds  encountered. 
Propagation  paths  that  are  longer  than  or  comparable  to  the 
size  of  the  clouds,  such  as  paths  1  and  2,  experience  regions 
of  weak  and  strong  turbulence  due  to  global,  and  to  a  lesser 
extent  due  to  intrinsic,  intermittency.  Paths  3-5  are  short 
compared  to  the  size  of  the  clouds  and  therefore  experience 
either  strong  or  weak  globally  intermittent  patches.  Since 
path  3  does  not  actually  pass  through  a  cloud,  sound  waves 
propagating  along  this  path  experience  very  little  scattering. 
Paths  4  and  5  are  through  similar  globally  strong  regions,  but 
scattering  is  stronger  along  path  4  because  of  differences  in 
intrinsic  intermittency. 

Several  previous  studies  have  dealt  with  the  effect  of 
intrinsic  intermittency  on  electromagnetic14'15  and 
acoustic 16-1 9  propagation.  The  underlying  idea  of  these  stud¬ 
ies  is  that  if  the  propagating  waves  sample  a  region  smaller 
than  the  outer  scale,  the  index-of-refraction  structure- 
function  parameter  C2  for  the  scattering  region  (which  de¬ 
scribes  the  strength  of  the  inertial-subrange  spectrum  in  this 
region)  becomes  a  random  variable.  The  statistics  of  the 
varying  C2  can  be  described  by  the  well-known  log-normal 
model.”  According  to  this  model,  the  smaller  the  sample 
region,  the  greater  the  fluctuations  in  C2. 

Other  studies  have  explicitly  considered  global  intermit¬ 
tency  effects  on  wave  propagation,  or  provide  treatments  for 
intermittency  effects  that  are  not  specifically  global  or 
intrinsic."  In  these  studies,  it  is  generally  assumed  that  the 
length  of  the  propagation  path  is  large  compared  to  the  outer 
scale.  Hence,  the  fluctuations  in  C2  resulting  from  intrinsic 
intermittency  are  unimportant.  However,  statistical  param¬ 
eters  for  the  turbulence,  including  C2,  may  still  vary  as  a 
result  of  globally  intermittent  processes.  Hentschel  and 
Procaccia"  considered  an  effect  of  intrinsic  intermittency, 
namely,  the  modification  of  the  well-known  -11/3  power 
law  for  inertial-subrange  turbulence,  over  global-scale  propa¬ 


•  27 

gation  paths.  Petenko  and  Shurygin"  combined  a  two- 
regime  model  for  global  intermittency  with  a  log-normal 
model  for  intrinsic  intermittency. 

In  comparison  to  previous  works,  in  this  article  we  study 
the  effects  of  combined  intrinsic  and  global  intermittency  on 
coherence  for  line-of-sight,  plane-wave  sound  propagation 
through  a  turbulent  atmosphere.  Of  course,  in  order  to  model 
the  propagation  in  intermittent  turbulence,  one  must  first 
model  the  intermittency.  A  sizable  portion  of  the  article  is 
concerned  with  a  quasiwavelet  (QW)  model  of  intermittency. 
A  QW  is  a  localized  perturbation,  analogous  to  a  turbulent 
eddy,  to  the  temperature  and/or  velocity  field."  ’  Several 
QW  models  of  atmospheric  turbulence  have  been  developed 
recently,  in  which  turbulence  is  represented  as  a  collection  of 
self-similar  QWs  of  many  different  sizes.  ’  These  and 
similar  models  have  been  used  for  studies  of  wave  propaga- 
tion  and  scattering  in  a  turbulent  atmosphere."  Due  to 
their  spatially  localized  nature,  QWs  can  be  a  useful  tool  for 
describing  turbulence  and  other  random  phenomena  with  in¬ 
trinsic  and  global  intermittency  features.  QWs  are  also  espe¬ 
cially  well  suited  to  modeling  inhomogeneous,  anisotropic 
turbulence  because,  unlike  customary  wavelets,  they  may  be 
distributed  and  oriented  nonuniformly  in  space  according  to 
any  desired  joint  probability  distributions. 

In  this  paper,  we  first  obtain  formulas  for  the  three- 
dimensional  (3D)  spectrum,  correlation  function,  variance, 
and  kurtosis  predicted  by  a  QW  model  of  temperature  fluc¬ 
tuations  with  intrinsic  intermittency.  Then,  the  QW  model  is 
generalized  to  account  for  global  intermittency.  Although 
both  temperature  and  velocity  fluctuations  affect  the  coher¬ 
ence  of  a  sound  wave  propagating  in  a  turbulent  atmosphere, 
their  contributions  are  additive  and  hence  can  reasonably  be 
studied  separately.6,28,30  We  therefore  simplify  the  develop¬ 
ment  of  a  QW  model  of  turbulence  with  intrinsic  and  global 
intermittency  by  considering  only  temperature  fluctuations. 
Reference  32  reports  on  preliminary  results  in  developing  a 
QW  model  of  intermittent  velocity  fluctuations.  Also,  al¬ 
though  we  apply  the  model  here  to  sound-wave  propagation 
through  turbulence,  it  potentially  has  broader  applications, 
e.g.,  to  electromagnetic  propagation,  geological  heterogene- 

38 

ities  close  to  Earth’s  surface  and  within  the  Earth, 
boundary-layer  meteorology,  and  eddy-based  simulation  of 
turbulence. 

This  paper  is  organized  as  follows.  In  Sec.  II,  a  QW 
model  of  temperature  fluctuations  with  intrinsic  intermit¬ 
tency  is  developed.  The  model  is  generalized  in  Sec.  Ill  to 
include  global  intermittency  along  a  one-dimensional  path. 
In  Sec.  IV,  this  statistical  framework  for  describing  intrinsic 
and  global  intermittency  is  used  as  the  basis  for  a  theory  of 
the  coherence  function  of  a  sound  wave  propagating  in  a 
turbulent  atmosphere.  Predictions  of  the  theory  are  studied 
with  numerical  simulation  methods,  and  analytical  approxi¬ 
mations  are  derived  and  compared  to  the  simulations.  In  Sec. 
V,  the  obtained  results  are  summarized. 

II.  QW  MODEL  FOR  INTERMITTENT  SCALAR 
FLUCTUATIONS 

In  this  section,  basic  formulas  describing  a  QW  model 
for  fluctuations  in  scalar  quantities  are  briefly  discussed. 
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(The  reader  may  refer  to  articles  cited  in  the  Introduction  for 
more  details.)  Then,  these  formulas  are  used  to  develop  QW 
models  of  temperature  fluctuations  with  intrinsic  and  global 
intermittency. 

A.  Basic  formulas 

The  underlying  idea  of  the  QW  model  is  to  represent  a 
random  field  as  a  collection  of  randomly  placed  and  oriented 
objects.  The  objects,  like  customary  wavelets,  are  spatially 
localized  and  based  on  rescaling  and  translation  of  a  parent 
function.  However,  some  properties  of  customary  wavelets, 
such  as  zero  mean,  can  be  relaxed.  Their  sizes  and  ampli¬ 
tudes  can  be  selected  according  to  any  desired  scaling  laws. 
Lovejoy  and  Mandelbrot  have  called  this  sort  of  self¬ 
similar  representation  a  “fractal  sum  of  pulses.”  Although 
most  of  the  discussion  to  follow  explicitly  deals  with  turbu¬ 
lent  temperature  fields,  the  same  modeling  approach  can  be 
applied  to  other  types  of  scalar  fluctuations.  For  a  QW  model 
of  turbulence,  the  individual  QWs  are  roughly  analogous  to 
turbulent  eddies. 

A  single  temperature  QW  is  a  perturbation  A'Tm'(R)  to 
the  mean  field  given  by 

A7""(R)  =  7a"A7’(/(|R  -  b“"|/a„) .  (1) 

Here,  R=(x,y,z)  are  the  Cartesian  coordinates,  and  a 
=  1,2, ...  ,N  and  n  =  l,2, ...  ,Na  are  two  indices,  which  de¬ 
termine  a  particular  QW.  The  index  a  determines  the  size  aa 
of  a  QW:  There  are  N  QW  sizes  arranged  in  a  diminishing 
order:  al> a2>  ■  ■  ■  > aN.  The  largest  a1  and  smallest  aN  sizes 
are  of  the  order  of  outer  and  inner  scales  of  turbulence.  In  a 
QW  model,  there  are  Na  QWs  with  the  same  size  an\  indi¬ 
vidual  QWs  in  a  size  class  are  distinguished  by  the  index  n. 
Furthermore,  in  Eq.  (1)  b“"  are  the  coordinates  of  the  center 
of  the  an th  QW,  A Ta  is  its  amplitude,  t“"  is  a  random  sign 
factor  with  zero  mean  and  unit  variance,40  and  /  is  the  parent 
function  describing  the  shape  of  the  QW.  Note  that  b“"  is  a 
random  vector  within  the  volume  V  where  turbulence  is 
modeled  and  /  is  the  same  for  all  QWs.  Different  parent 
functions  /  can  be  used  in  QW  models.  In  what  follows,  a 
theoretical  development  will  be  done  for  an  arbitrary  parent 
function  /.  Specific  results  will  subsequently  be  obtained  for 
the  Gaussian  parent  function,  given  by  f(x) 
=  (2tt)312  exp(-;^/2),  where  y  is  a  nondimensional  argu¬ 
ment.  The  overall  field  of  temperature  fluctuations  T{ R)  is 
the  sum  of  contributions  from  the  individual  QWs  as  fol¬ 
lows: 

N  Na 

T(R)=  2  S  A7’“”(R).  (2) 

a=  1  n=\ 

Here,  the  sums  are  taken  over  all  QW  sizes  used  and  over  all 
QWs  within  a  particular  size  class. 

A  correlation  function  B{ R)  of  temperature  fluctuations 
is  defined  as  follows: 

B(R,-R2)  =  (f(RI)f(R2)),  (3) 

where  the  brackets  ( )  denote  ensemble  average  and  it  is 
assumed  that  temperature  fluctuations  are  statistically  homo- 
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geneous.  The  3D  spectrum  <F( k)  of  temperature  fluctuations, 
where  k  is  the  turbulence  wave  vector,  is  a  Fourier  transform 
of  B( R).  For  the  case  of  isotropic,  homogeneous  turbulence, 
<!>(  k)  depends  only  on  the  magnitude  of  the  vector  k.  Sub¬ 
stituting  Eqs.  (1)  and  (2)  into  Eq.  (3),  and  assuming  that  the 
b“"  and  r“"  are  mutually  independent,  the  following  formula 
for  <!>(/<•)  can  be  obtained  by  averaging  over  an  ensemble  of 
random  realizations: 30 

N 

O(ac)  =  8 77 3 2  ^^A7^F2(/Cflj.  (4) 

a=  1 

Here,  (f>a=N cfc’J'V  is  the  packing  fraction  of  the  QWs  with 
the  same  size  aa,  which  is  proportional  to  the  ratio  of  the 
volume  occupied  by  these  QWs  to  the  total  volume  V  where 
temperature  fluctuations  are  simulated.  Furthermore,  in  Eq. 
(4)  F(t;)  is  the  spectral  parent  function;  that  is,  the  Fourier 
transform  of  the  parent  function  fix)-  For  the  Gaussian  par¬ 
ent  function,  T(Q  =  exp(-£2/2). 

The  correlation  function  B(R)  is  the  inverse  Fourier 
transform  of  .  By  the  convolution  theorem,  the  inverse 
transform  of  F2(£)  is  the  convolution  of  f(x)  with  itself. 
Hence,  the  correlation  function  is 

N 

B(R)  =  'E  ct>^T2J2(Rlaa),  (5) 

a=\ 

where  /2(;^)  is  defined  as 

flix)  =/(*)  *f(x)  =  J  d3x'f(\x-x'\)f(x')-  (6) 

The  preceding  formulas  are  quite  general  and  can  be 
used  to  model  random  fields  having  a  range  of  properties.  To 
capture  the  properties  of  turbulence,  certain  scaling  relation¬ 
ships  should  be  imposed  on  the  QW  sizes  aa,  the  packing 
fractions  cf>a ,  and  the  temperature  amplitudes  A  Ta.  For  ex¬ 
ample,  to  model  statistically  isotropic,  homogeneous,  and 
nonintermittent  temperature  fluctuations,  the  scaling  relation¬ 
ships  can  be  chosen  as  follows:30 

aa  =  4>a  =  const,  AT„=  A^  {aja^) 1/3 , 

(7) 

where  /a  is  a  positive  parameter  usually  much  smaller  than  1 . 
According  to  the  first  of  these  scaling  relationships,  the 
neighboring  sizes  have  a  constant  ratio  (e_At),  which  mimics 
a  self-similar  turbulent  cascade  of  eddies  with  different  sizes. 
The  second  relationship  in  Eq.  (7)  ensures  that  QWs  with 
different  sizes  occupy  the  same  total  volume,  which  is  ex¬ 
pected  for  nonintermittent  turbulence.  Finally,  the  third  rela¬ 
tionship  produces  the  classical  Kolmogorov  spectrum  of 
temperature  fluctuations  in  the  inertial  subrange  as  should  be 
the  case  for  isotropic,  homogeneous  turbulence.  After  substi¬ 
tution  of  Eq.  (7)  into  Eq.  (4)  and  some  algebra,  a  spectrum 
<D(/r)  qualitatively  similar  to  von  Karman’s  is  obtained. 30 

Figure  2(a)  shows  an  example,  random  QW  field.  This 
realization  was  made  from  362  size  classes  ranging  from 
aN=  0.2  m  to  «]  =  10  m.  (QWs  smaller  than  0.2  m  would  not 
be  apparent  in  the  image.)  The  packing  fraction  <f>  is  0.0023 
for  all  size  classes  and  AT]  =0.5  K.  The  QWs  are  distributed 
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FIG.  2.  Examples  of  synthetic  turbulence  fields  produced  by  random  QW 
ensembles,  (a)  Example  without  intermittency.  (b)  Example  with  intrinsic 
intermittency  (A. =0.5).  (c)  Same  as  (b),  except  that  the  packing  fraction  has 
been  multiplied  by  1/3,  and  the  amplitudes  by  y3.  All  scales  are  in  K. 


uniformly  in  a  160X60X80  m3  volume.  A  100X20  m2 
cross  section  along  the  xz-plane  is  shown.  To  mitigate  edge 
effects,  the  cross  section  is  centered  within  the  volume  such 
that  all  displayed  points  are  at  least  30  m  (3«|)  from  a  side 
of  the  volume. 

6.  Intrinsic  intermittency 

1.  Scaling  relationships  and  cascades 

To  model  temperature  fluctuations  with  intrinsic  inter¬ 
mittency,  the  scaling  relationships  for  cia,  <pa,  and  A Ta  must 
be  chosen  to  capture  the  main  properties  of  such  fluctuations. 
Regardless  of  whether  intermittency  is  present,  we  anticipate 
that  the  turbulence  is  a  self-similar  cascade  process  and  thus 
the  distribution  of  size  classes  continues  to  be  dictated  by  the 
first  relationship  in  Eq.  (7).  However,  for  turbulence  with 
intrinsic  intermittency,  small  eddies  become  less  and  less 
space  filling1”  so  that  the  second  relationship  in  Eq.  (7)  does 
not  hold.  To  accommodate  this  property  of  intermittent  scalar 
turbulence,  we  assume  (similarly  to  Ref.  32  where  intermit¬ 
tent  velocity  fluctuations  were  considered)  that  the  packing 
fraction  4>a  decreases  as  the  QW  size  becomes  smaller  as 
follows: 

4>a=  <KaJaxY.  (8) 

Here,  (f>  is  a  constant  that  is  equal  to  the  packing  fraction  (/>, 
of  the  largest  QWs  and  A.S=0  is  a  parameter  characterizing  a 
degree  of  intrinsic  intermittency.  The  greater  the  parameter 
X,  the  more  intermittent  the  turbulence;  A  =  0  corresponds  to 
the  nonintermittent  case. 

The  cascade  mechanism  by  which  the  small  eddies  be¬ 
come  less  space  filling  controls  the  spatial  distribution  of  the 
QWs.  In  turbulence  literature,  the  cascade  mechanism  is  of¬ 
ten  conceived  as  a  process  repeated  from  one  “generation”  to 
the  next:  Parent  eddies  produce  children  eddies  that  are  con¬ 
tained  within,  yet  do  not  entirely  fill,  the  space  of  the  parents. 
For  example,  a  parent  eddy  of  size  a  may  break  down  into 
several  eddies  of  size  al  2.  In  three  dimensions,  2 3  children 
eddies  would  be  needed  to  maintain  the  packing  fraction. 
Less  than  23  eddies  would  decrease  the  packing  fraction. 


This  description  coincides  to  the  well-known  beta  model  of 
Frisch  et  al.,  ~  which  defines  the  parameter  f3  as  the  number 
of  children  eddies  divided  by  23.  Since  [i=  <ba+\  /  </>„ 
=  {aa+ll aa)x  when  aa+l/aa=l/2,  it  follows  that  f3=  2~\ 
Frisch  et  al . 12  showed  that  (3=  2°~3,  where  D~  2.5  is  the 
fractal  dimension.  Hence,  \  =  3-Z)~0.5. 

Schmitt  et  al.4]  described  a  procedure  for  “scale  densi- 
fication”  of  the  beta  model  that  removes  the  restriction  of  the 
parent  eddy  size  being  an  integer  multiple  of  the  size  of  the 
children.  Lovejoy  and  Scherzer  generalized  the  beta  model 
so  that  regions  are  not  perfectly  active  or  inactive.4"  Con¬ 
ceivably,  the  beta  model  or  such  a  generalization  thereof 
could  be  adapted  to  place  the  QWs  in  space  according  to  a 
cascade  process.  However,  there  are  many  challenges,  such 
as  how  to  provide  a  continuous  placement  of  the  QWs  and 
how  to  derive  statistical  results  when  the  size  classes  are  not 
independent.  Also,  the  beta  model  is  highly  idealized:  Actual 
eddies  are  swept  through  space  during  a  cascade  process  and 
the  cascade  process  may  be  at  different  ages  in  different 
regions  of  space.  For  present  purposes,  we  do  not  attempt  to 
formulate  a  cascade  spatial  construction  with  the  QWs;  the 
positions  of  QWs  in  each  size  class  are  independent  of  each 
other  and  of  the  positions  of  QWs  in  the  preceding  size  class. 

This  assumption  is  actually  consistent  with  the  fractal  sum  of 

39 

pulses  method  described  by  Lovejoy  and  Mandelbrot/ 

Intrinsic  intermittency  of  turbulence  can  also  affect  the 
scaling  relationship  for  QW  amplitudes  A Ta.  To  accommo¬ 
date  this  possibility,  the  scaling  relationship  for  A Ta  is  speci¬ 
fied  as  follows: 

ATa=  ATMJai)V\<f>l4>a)v=  (9) 

where  77  is  another  parameter  characterizing  the  intrinsic  in¬ 
termittency.  If  77=  0,  the  scaling  relationship  described  by  Eq. 
(9)  coincides  with  Eq.  (7).  The  value  of  the  parameter  77 
depends  on  the  problem  under  consideration.  For  example, 
the  value  of  77  for  turbulent  temperature  fluctuations  might 
differ  from  that  for  random  heterogeneities  within  Earth.  For 
turbulence,  a  value  for  77  can  be  derived  from  conservation 
considerations,  as  has  been  previously  demonstrated  for  a 
QW  model  of  turbulent  velocity  fluctuations  with  intrinsic 
intermittency.  "  Here,  we  extend  this  treatment  to  conserva¬ 
tive  scalar  fluctuations  such  as  temperature. 

First,  consider  the  specific  (per  unit  mass)  turbulent  ki¬ 
netic  energy  (TKE).  The  TKE  for  size  class  a  is  proportional 
to  (f>av2a,  where  va  is  the  characteristic  velocity  scale.  The 
rate  of  transfer  of  TKE  from  eddies  with  scale  aa  to  smaller 
scales  should  be  proportional  to  the  TKE  divided  by  the  time 
scale  for  the  size  class,  ajv  a.  Since  kinetic  energy  is  neither 
created  nor  destroyed  within  the  inertial  subrange,  the  rate  of 
transfer  should  be  invariant  for  turbulence  in  equilibrium. 
Hence,  fjx’J aa=  <pv]l al  for  all  a.  Substituting  with  Eq.  (8), 
we  then  have 

va=vl(ajal){l~x)l3.  (10) 

Similarly,  the  temperature  variance  associated  with  size  class 
a  is  <paA T2a.  If  the  rate  of  variance  transfer  is  to  be  preserved 
(as  it  would  be  for  a  conservative  scalar  in  steady  turbu¬ 
lence),  we  must  have  <f>aAT2av al aa=  4>AT\v \l ax.  Substituting 
with  Eqs.  (8)  and  (9)  then  yields 
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Ara=A  r1(fla/fl1)(1-x,/3  (li) 

and  hence  77  in  Eq.  (9)  is  1/3.  The  value  ?7>0  implies  that 
the  eddies  (QWs)  must  spin  relatively  faster,  and  have  a 
stronger  temperature  amplitude,  to  compensate  for  the  de¬ 
creasing  packing  fraction. 

Figures  2(b)  and  2(c)  are  similar  to  Fig.  2(a),  except  that 
intrinsic  intermittency  with  A.  =  0.5  and  77=  1/3  is  included. 
Figure  2(c)  furthermore  decreases  the  packing  fraction  </>  by 
1/3  (to  0.000  76)  for  all  size  classes  while  increasing  A  7]  by 
\  3  (to  0.866  K);  according  to  Eq.  (4),  these  rescalings  do  not 
alter  the  spectrum  or  correlation.  With  intrinsic  intermittency, 
there  are  noticeably  fewer  small  QWs,  although  they  are 
stronger.  Decreasing  packing  fraction  also  enhances  the  ap- 
pearence  of  intermittency. 


greater  than  the  inner  scale  of  turbulence,  which  is  of  order 

aN- 

Although  Eq.  (16)  is  specific  to  the  Gaussian  parent 
function,  results  for  other  reasonable  parent  functions  have  a 
similar  dependence  on  wave  number."  ’  The  scaling  laws 
[Eq.  (7)  for  classical  turbulence,  or  Eqs.  (8)  and  (9)  for  in¬ 
termittent  turbulence],  when  combined  with  the  fractal  size 
distribution,  lead  to  inertial  subranges  with  a  slope  indepen¬ 
dent  of  the  parent  function.  The  choice  of  parent  function  is 
mainly  significant  in  the  energy-containing  subrange  ( kcix 
<8  1)  and  in  the  transition  between  subranges. 

As  for  the  correlation  function,  one  can  show  for  the 
Gaussian  parent  function  that 

/2(A')  =  87T9/2e-^/4.  (17) 


2.  Spectra  and  correlations 

With  these  revised  scaling  relationships,  we  can  now 
calculate  the  spectrum  Tf k)  from  Eq.  (4).  In  this  equation, 
cf>a  and  A Ta  are  replaced  with  their  values  given  by  Eqs.  (8) 
and  (9),  respectively,  and  aa  is  replaced  with  its  value  given 
by  the  first  relationship  in  Eq.  (7).  In  the  resulting  formula, 
assuming  that  /r<Sl,  the  sum  over  a  is  replaced  with  the 
integral  over  a  using  Eq.  (8)  from  Ref.  30.  As  a  result,  we 
obtain  a  formula  for  the  3D  spectrum  of  temperature  fluctua¬ 
tions  with  intrinsic  intermittency  as  follows: 

rka  1 

<s>{K)=c{Kai)-im-v\  ei3+vF2im.  (12) 

J  kaN 


Here,  the  new  parameter  v  is  a  combination  of  the  intrinsic 
intermittency  parameters  A.  and  77, 


17=  A(1  —  2 77)  ~  1/6, 
and  the  coefficient  C  is  given  by 


C  =  8tt 


,  1  at2 

P 


(13) 


(14) 


Equation  (12)  describes  a  3D  spectrum  of  temperature  fluc¬ 
tuations  with  intrinsic  intermittency  for  an  arbitrary  parent 
function  F(£).  If  v=0,  this  spectrum  coincides  with  that  for 
nonintermittent  temperature  fluctuations  [Eq.  (9)  in  Ref.  30]. 
Similarly,  from  Eq.  (5)  one  finds 


B(R)  = 


4>m\!  r  \2/3+"  fWaN 

P  \aj  JR/ai 


X~5li~vf2(x)dx- 


(15) 


For  the  Gaussian  spectral  parent  function  /*’(£)  =  ex p(-f2/2), 
the  integral  on  the  right-hand  side  of  Eq.  (12)  can  then  be 
calculated  as 


=  ~(Kai)  11/3  "[7(11/6+  v/2,  ica]) 

-  7(1 1/6  +  17/2,  o^)].  (16) 

Here,  7 (a,x)  is  the  incomplete  gamma  function  and  the  sub¬ 
script  G  stands  for  “Gaussian.”  For  the  case  of  nonintermit¬ 
tent  turbulence,  when  v=0,  Eq.  (16)  coincides  with  Eq.  (11) 
from  Ref.  30  as  it  should.  For  sound  propagation  in  a  turbu¬ 
lent  atmosphere,  the  sound  wavelength  is  nearly  always 


One  then  finds,  from  Eq.  (15), 
i?aC(  R  \ 2/3+11 

BC(R)  =  ~ ^  —  [IX-  1/3  -  vl2,R2/4af) 

2 a i  \ j 

-T(-  1/3-  v/2,R2/4a2N)],  (18) 


where  T (a , x)  is  the  complimentary  incomplete  gamma  func¬ 
tion.  Since  -1/3-  v/2  is  negative  for  v>-2/3,  and  the  com¬ 
plete  gamma  function  T(a)  is  undefined  for  a  negative  argu¬ 
ment,  it  is  helpful  to  apply  the  recursion  formula  T(a+ 1  ,x) 
=  aT(a,x)+xae~x  to  rewrite  Eq.  (18)  as 


BC(R) 


7 r^C 
(2/3  +  v)a\ 


e 


-R2!Aa\ 


R  y 2/3+11 

2a!/ 


XT(2/3  -  v/2,R2/4a2l)  - 


aN 

a1 


2/3+1' 

£-R~/4a 


2 

N 


+ 


2/3+v 

T(2/3  -  vl2,R2l4a2N) 


(19) 


Temperature  fluctuations  with  scale  less  than  aN  often  do 
not  affect  the  coherence  function  and  other  statistical  mo¬ 
ments  for  line-of-sight  sound  propagation.9  Therefore,  in 
Eqs.  (16)  and  (19),  it  is  often  reasonable  to  set  aN=  0.  The 
terms  involving  aN  thus  vanish.  In  the  energy  subrange  of 
turbulence,  where  Kal<l,  the  incomplete  gamma  function 
in  Eq.  (16)  can  be  approximated  as  follows:  7(1 1/6 
+  W2,A-2a2)~(A:ai)11/3+17(ll/6+i7/2).  In  this  subrange,  the 
3D  spectrum  of  temperature  fluctuations  does  not  depend  on 
k  and  is  given  by  <I>G=C/(1  1/3  +  v).  In  the  inertial  subrange, 
where  Kci\i>  1,  the  incomplete  gamma  function  in  Eq.  (16) 
can  be  replaced  with  its  asymptotic  value  for  large  values  of 
the  argument.  As  a  result  (with  aN=  0),  we  have 


CT(  11/6 +17/2),  x_im_ 
- - - (kci  j) 


(20) 


It  follows  from  this  formula  that  increasing  intermittency 
(increasing  v)  steepens  the  decay  of  ^(/c)  in  the  inertial 
subrange.  It  is  worthwhile  to  compare  <f>c(/r)  given  by  Eq. 
(20)  with  the  Kolmogorov  spectrum 
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(21) 


<b(K)=AC2rK-'m, 


which  is  valid  in  the  inertial  subrange  of  nonintermittent  tur¬ 
bulence.  Here,  A  — 0.033  is  a  numerical  constant  and  C\  is 
the  structure-function  parameter  for  temperature  fluctuations. 
For  v=0,  Eq.  (20)  has  the  same  /e-dependence  as  the  classi¬ 
cal  (nonintermittent)  Kolmogorov  spectrum.  If  we  constrain 
the  parameters  in  the  QW  representation  such  that 


Ca 


-11/3 

1 


,3^  2A 

%af  “r(l  1/6+ v/2)Cr’ 


(22) 


Eq.  (22)  implies  A 7" x  —  a J/3.  Equation  (20)  now  becomes 


<J>G(/c)  =ACjK  U  i(kal)  ". 


(23) 


This  result  is  consistent  with  Eq.  2.20  of  Hentschel  and 
Procaccia.  “  In  the  inertial  subrange,  the  intermittency  effect 
modifies  the  /c-ll/3  power  law  to  k~1  1/3~"~  k~2316. 

For  Rlax^  1,  the  complementary  incomplete  gamma 
function  in  Eq.  (18)  can  be  replaced  with  the  complete 
gamma  function,  and  one  has  for  the  correlation  function 
(neglecting  the  effect  of  aN) 


BC(R ) 


tt212 C 
(2/3  +  v)a3{ 


(- 

\2al 


2/3+i* 

T(2/3  -  v/2) 


(24) 


The  structure  function  for  the  field  is  defined  as  D(R) 
=  2[B(0)  -  B(R)].  For  small  separations,  the  structure  function 
becomes 


773/2Cr(2/3- v/2)/  R  Y+213 
C  (1/3  +  v!2)a\  \2aJ 

Substituting  with  Eq.  (22),  we  have 
DC(R)  «  C2^2/3(^)  * 
where  we  have  made  the  definition 


(25) 


(26) 


(1/3  +  iV2)r(ll/6  +  v!2) 
2i/3-V7J3/2  y(2I3-VI2) 


(27) 


(With  this  definition.  A  — 0.033  when  v=Q,  as  expected.)  If 
aN  is  not  identically  zero,  we  have  for  R<aN 


TTinC  i  R  \"+2/3 

BC(R)  =  Ty  —  [r(-  1/3  -  vl2,R2IAa\ ) 

2a j  \2ax! 

-T(-  1/3-  v/2,R2IAa2N)].  (28) 


Using  Eq.  (20),  the  3D  spectrum  <f>G(/c)  [normalized  by 
C/(ll/3),  the  value  of  <FC(0)  for  n=0]  is  plotted  in  Fig.  3 
versus  the  normalized  wave  parameter  ku  ,  for  v=0,  0.25, 
and  0.5.  The  spectrum  <J>G(/c)  has  two  distinct  regions.  It  is 
almost  constant  in  the  energy  subrange,  where  kci\<  1,  and 
has  a  power  law  dependence  on  k  in  the  inertial  subrange, 
where  Kaxi>  1.  This  is  consistent  with  asymptotic  behavior 
of  the  spectrum  in  the  energy  and  inertial  subranges  consid¬ 
ered  above.  Furthermore,  it  follows  from  Fig.  3  that  the 
greater  the  parameter  v ,  the  smaller  are  the  values  of  the 


FIG.  3.  (Color  online)  Normalized  spectrum  <I>G  of  intermittent  temperature 
fluctuations  vs  the  normalized  wave  number  kci  ]  for  different  values  of  the 
intrinsic  intermittency  parameter  v. 


spectrum  ^(/c).  This  tendency  is  more  pronounced  in  the 
inertial  subrange  where  the  spectral  slope  steepens  with  in¬ 
creasing  value  of  this  parameter. 

Figure  4  compares  theoretical  predictions  to  the  struc¬ 
ture  function  estimated  from  random  realizations  of  QW 
fields.  The  case  considered  is  the  same  as  Fig.  2(b),  namely, 
aN=  0.2,  a]  =  10m,  </>=  0.0023,  and  AT]  =0.5  K.  The  esti¬ 
mates  were  derived  from  correlation  functions  of  256  ran¬ 
dom  realizations.  The  theoretical  curves  shown  are  the  exact 
result  (with  discrete  size  classes)  based  on  Eq.  (5),  the  con¬ 
tinuous  size-class  approximation  based  on  Eq.  (18),  and  the 
inertial-subrange  approximation,  Eq.  (25).  The  simulations 
and  exact  result  are  nearly  indistinguishable.  For  very  small 
separations,  the  continuous  approximation  deviates  some¬ 
what  from  the  exact  result.  The  inertial-subrange  approxima- 


FIG.  4.  Structure  function  for  turbulence  with  intrinsic  intermittency.  The 
solid  line  is  an  exact  curve  with  discrete  spacing  between  size  classes;  the 
dashed  line  is  the  corresponding  approximation  for  continuous  size  classes. 
The  dotted  line  is  an  inertial  subrange  approximation.  The  circles  show  the 
actual  structure  functions  determined  from  256  random  QW  realizations. 
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tion  is  seen  to  be  poor  for  this  particular  case;  apparently,  aN 
is  insufficiently  small  compared  to  a1  for  this  approximation 
to  be  accurate. 


3.  Variance  and  kurtosis 

The  variance  of  temperature  fluctuations  cr  is  needed  in 
many  applications.  For  both  intermittent  and  nonintermittent 
turbulences,  cr  can  be  determined  by  integrating  the  3D 
spectrum  of  temperature  fluctuations  or  by  taking  the  limit 
/?— > 0  of  Eq.  (15).  Using  the  latter  approach,  we  have  from 
Eq.  (5) 

O2  =  2  Oa,  (29) 

a- 1 


where  the  contribution  to  the  variance  from  size  class  a  is 

oi=cf>aAT2J2(0)-  (30) 


To  obtain  the  variance  for  the  intermittent  case,  4>a  and  A  Ta 
in  Eq.  (29)  are  replaced  with  their  values  given  by  Eqs.  (8) 
and  (9),  respectively.  We  thus  obtain 

a2a=cf>AT\{acJal)m+v.  (31) 


Setting  aa=axe  l±<a  1 1  and  recognizing  Eq.  (29)  now  as  a 
geometric  series,  we  explicitly  determine  the  sum  as 


cr  =  ^.Ar^/2(0) 


1  _  g-M2/3 +v)N 
1  _  e-/i(2/3+i>) 


0A7fc(O) 

/r(  2/3  +  v) 


(32) 


The  second  approximate  form  is  valid  when  there  are  many 
closely  spaced  size  classes  (/rA8>  1  and  /r<S  1).  Applying 
Eqs.  (14)  and  (17)  yields  the  variance  for  Gaussian  QWs  as 
follows: 


2  tt2'2  C 

G~  (2/3  +  v)a\' 


(33) 


We  see  that  the  variance  decreases  with  increasing  intermit- 
tency.  When  v=  1/6,  cr  is  4/5  times  its  value  without  inter- 
mittency. 

The  kurtosis  K  is  defined  as  the  normalized  fourth  mo¬ 
ment  of  a  random  field,  namely,  K=(T4)I cr4.  Since  intermit- 
tency  generally  leads  to  a  kurtosis  larger  than  the  value  for  a 
Gaussian  random  variable,  namely,  K=  3,  determination  of 
the  kurtosis  is  of  much  interest  (e.g.,  Refs.  43  and  44).  A 
general  formula  for  (T4)  in  the  QW  model  is45 

N  j  N  \2 

(T4)  =  Kt'£  +  3  2  a2  \  ,  (34) 

a— 1  \  a— 1  / 


where  fl  =  and  K,  is  the  fourth  moment  of  the  r“" 

(which  as  described  earlier,  have  a  unit  variance).  For  a  QW 
model  in  which  the  r™  are  ±1  with  equal  probability,  K, 
=  1 .  For  a  model  in  which  the  t“"  are  normally  distributed, 
K,= 3.  Equation  (34)  is  valid  for  both  intermittent  and  non¬ 
intermittent  turbulence.  To  obtain  the  value  of  K  for  the  in¬ 
termittent  case,  aa,  cf>a,  and  A7  „  are  replaced  with  their  val¬ 
ues  given  by  Eqs.  (7)-(9),  respectively.  We  thus  obtain 


FIG.  5.  Coefficient  6,  which  controls  the  kurtosis  of  the  QW  field,  vs  the 
intrinsic  intermittency  parameters  X.  and  77. 


cfjATlllK.  " 

K=  3  +  - - \ - -2j  e" 


(a-\)q 


(35) 


where  q  =  \(\  -Atj).  Unlike  the  variance,  the  fourth  moment 
[Eq.  (34)]  is  not  a  simple  sum  of  contributions  from  the  size 
classes.  However,  Eq.  (35)  shows  that  the  deviation  from 
Gaussian  statistics,  K—  3,  can  be  considered  as  the  sum  of 
contributions  from  the  size  classes  if  the  variance  is  fixed.  By 
calculating  the  sum  in  Eq.  (35),  a  formula  for  the  kurtosis  of 
temperature  fluctuations  with  intrinsic  intermittency  results 

<t)AT\VLK,  1  -  e-/44/3+9)v 

K=3+  ~A  !  _  e-/44/3+9) 

u(2/3  +  v)2VLK.  ,  , 

~3  +  tm  +  qXfMT  36 


The  second,  approximate  form  applies  when  fiN  §>  1  and  /i 

The  preceding  formulas  are  valid  for  any  parent  func¬ 
tion.  For  the  Gaussian  parent  function,30  n  =  (7r/2)3/2(27r).6 
Combining  this  result  with  Eq.  (17),  we  obtain  the  kurtosis 
for  the  case  of  the  Gaussian  parent  function 


K  o  ,  ^>K' 

G  3(2t t)3/2</>' 


(37) 


Here,  the  coefficient  9  is  a  combination  of  parameters  v  and 
q  (and  thus  A.  and  rj)  as  follows: 

(1  +  3v/2)2  ,  x 

0=  - - -.  (38) 

1  +  3  g/4 


For  nonintermittent  turbulence,  when  A  =  0  and  0=  1,  Eq.  (37) 
coincides  with  Eq.  (44)  from  Ref.  30.  The  kurtosis  can  be 
enhanced  either  by  a  low  density  of  QWs  (small  </>//+)  or 
intrinsic  intermittency  ( 0>  1 ).  Illustrations  of  the  former 
mechanism  can  be  found  in  Ref.  30. 

In  Fig.  5,  the  coefficient  0  is  plotted  versus  the  param¬ 
eters  A  and  rj.  It  follows  from  this  figure  that  8  varies  in  the 
range  from  1  to  2.23  for0^A^0.5  and  0«  975=  0.5.  Figure  5 
can  be  used  in  choosing  parameters  of  a  QW  model  of  tem¬ 
perature  fluctuations  leading  to  a  particular  value  of  Kc. 
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To  test  the  preceding  expressions,  we  compare  them  to 
variance  and  kurtosis  from  simulations.  The  cases  considered 
are  the  same  as  shown  in  Fig.  2;  results  from  256  such  ran¬ 
dom  realizations  were  averaged.  For  the  case  without  intrin¬ 
sic  intermittency  and  0=  0.0023,  as  illustrated  in  Fig.  2(a), 
the  simulation  results  (with  theoretical  values  in  parentheses) 
were  (4=0.420  (0.420)  and  KG=?>A\  (3.34).  For  the  case 
with  intrinsic  intermittency  and  0=0.0023,  as  illustrated  in 
Fig.  2(b),  the  simulation  results  were  crc= 0.343  (0.350)  and 
Kc= 3.66  (3.56).  For  the  case  with  intrinsic  intermittency  and 
0=0.000  76,  as  illustrated  in  Fig.  2(c),  the  simulation  results 
were  (4=0.350  (0.350)  and  KG=4.S2  (4.68). 


C.  Global  intermittency 

As  discussed  earlier,  global  intermittency  involves  exter¬ 
nally  imposed  variations  in  turbulence  intensity  over  regions 
larger  than  the  outer  scale  of  the  turbulence.  In  this  section, 
we  consider  an  approach  to  incorporating  global  intermit¬ 
tency  into  a  QW  model. 

The  volume  V  where  the  turbulence  is  modeled  is  con¬ 
ceptually  subdivided  into  subvolumes  V,  with  characteristic 
scales  larger  than  a , .  Within  each  of  these  subvolumes,  QWs 
occur  as  described  in  Sec.  II  B  1,  although  the  parameters 
may  vary  from  one  subvolume  to  another.  In  principle,  any 
or  all  of  the  parameters  A Tb  aG  aN,  cf>,  p,  \,  and  77  could 
vary.  Variations  in  the  parameters  A  Tl  and  a  t  are  of  particu¬ 
lar  interest  since  they  depend  on  the  external  mechanism 
creating  the  turbulence.  The  inner  scale  aN  may  vary  but 
does  not  strongly  affect  sound  waves.  The  ratio  0/  / jl  and  the 
parameter  X  are  expected  to  have  a  physical  meaning  (the 
packing  of  QWs  as  a  function  of  their  size)  but  be  fixed  for  a 
particular  cascade  process,  such  as  turbulence.  (The  indi¬ 
vidual  parameters  0  and  p  are  constructive  in  the  QW  model 
but  should  not  vary  independently  for  a  particular  cascade 
process.)  According  to  Eq.  (14),  variations  in  AT)  and  a , 
result  in  random  values  of  the  coefficient  C.  Then,  it  follows 
from  Eqs.  (22)  and  (33)  that  the  values  of  C2r  and  (0.  are  also 
random  and  proportional  to  each  other. 

To  proceed,  we  need  statistical  models  for  the  variations 
in  AT,  and/or  ah  or  for  quantities  dependent  on  them,  such 
as  C.  Many  and  varied  statistical  models  for  global  intermit¬ 
tency  have  been  considered  in  literature.  Mahrt,11  Petenko 
and  Shurygin,  and  other  authors  used  dichotomous  regions 
of  strong  and  weak  activity.  Antonia  et  a/.46  considered  a 
linear  ramp  model.  Tatarskii  and  Zavorotnyi  considered  a 
gamma  probability  density  function  (pdf),  and  Frehlich2 1  a 
log-normal  pdf,  for  C2.  This  diversity  of  approaches  to  mod¬ 
eling  global  intermittency  reflects  the  different  mechanisms 
producing  it  as  well  as  the  challenges  of  describing  the  un¬ 
derlying  physics  with  tractable  models. 

In  the  following,  we  consider  a  basic  model  for  global 
intermittency  that  can  be  readily  related  to  a  QW  construc¬ 
tion.  Depending  on  how  the  parameters  are  adjusted,  the 
model  can  be  made  similar  to  most  previous  treatments  of 
global  intermittency.  It  involves  allowing  the  quantity  C  to 
vary  with  the  coordinate  x  along  a  linear  path;  in  our  case, 
this  path  is  the  propagation  path.  The  random  field  C(x)  is 
written  as 


C(x)  =  (C)  +  C'(x)  =  <C>[  1 +  £(*)],  (39) 

where  (C)  is  the  mean  value  and  C'(x)  is  its  fluctuating  part, 
and  £(x)  =  C'(x)/(C). 

Initially,  we  consider  a  two-state  Markov  process  for 
C(x).  (The  reader  may  refer  to  a  text  such  as  Wilks4  for  an 
introduction  to  two-state  Markov  processes.)  Within  each 
subvolume  V,  along  the  propagation  path,  the  turbulence  is 
either  active  (present)  or  inactive  (absent).  The  inactive  state 
is  characterized  by  A7j=0  and  C= 0;  the  value  for  a ,  is  thus 
immaterial.  The  values  of  A70  C,  and  a  t  are  the  same  for  all 
active  subvolumes.  Let  us  designate  Ca  as  the  value  of  C 
when  in  the  active  state,  and  the  actual  value  of  C  within  V, 
as  Cj=CaXh  where  A-,-  is  a  random  variable  equal  to  0  when 
Vj  is  an  inactive  state  and  1  when  it  is  active.  The  probability 
of  transitioning  from  an  inactive  state  in  the  volume  V,  to  an 
active  state  in  Vi+l,  while  moving  along  the  propagation  path 
from  one  subvolume  to  the  next,  is  designated  as  pol.  The 
probability  of  remaining  in  an  active  state  is  designated  pn, 
and  so  forth.  Since  there  are  only  two  possible  outcomes  for 
a  transition  from  a  given  initial  state,  we  must  have  p00 
+p01  =  1  and  Pio+Pn=l.  This  means  that  only  two  of  the 
transition  probabilities  are  independent.  By  convention,  these 
are  usually  taken  to  be  p01  and  pu.  If  the  Markov  process  has 
positive  memory,  the  probability  of  transitioning  to  a  particu¬ 
lar  state  is  greater  if  the  process  is  already  in  that  state  (p00 
>p  10,  pn>Poi)-  Also  of  interest  are  the  unconditional  prob¬ 
abilities  of  occurrence  for  the  inactive  and  active  states,  des¬ 
ignated  7T()  and  77]  =  1  —  77q.  Since  77]  equals  the  sum  of  p0]770 
and  pi  |77|,  we  have 


Pen 

77]  =  - . 

1  +P01~PU 


(40) 


The  following  relationships  can  also  be  proven: 


(Xi)  =  {X2)=  77„ 

(41) 

<A;2)  =  <(A,-<A;))2)=  77077i, 

(42) 

( Xi+jX'i  )  =  77077,  (p  1]  Poi)J  • 

(43) 

Statistics  involving  C,  or  C[  follow  after  appropriate  scaling 
by  Ca,  e.g.,  (Cj)=CaTTl  and  (C\2)  =  C2aTTQirx.  Statistics  for 
follow  from  those  for  X-  after  scaling  by  tt]1.  Designating 
the  distance  between  the  subvolumes  as  t,  we  can  write  Eq. 
(43)  as 

(X'(x  +  Ax)X'(x))  =  77077] (pn  -p01)|A'|,<  =  770771e“|A-t|/L, 

(44) 

where  L=f71n[l/(pu-p0i)].  The  transition  probabilities  can 
be  determined  from  the  mean  activity  level  77]  and  the  cor¬ 
relation  length  by  L  as  follows: 

Pen  =  77] (1  -  e~(IL),  pu  =  77,  +  TT0e~tlL.  (45) 

Next  let  us  suppose  C(x)  is  the  average  of  M  indepen¬ 
dent,  identical  two-state  Markov  processes,  namely. 
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r  M 

C(x)  =  -J2x,(m), 

(46) 

where  each  of  the  A4”'1  follow  Eqs. 

(40)-(43).  We  then  find 

<C(x))  =  CflT71, 

(47) 

flc(Ax)-<C'(x+Ax)C'(x)>- 

(48) 

B(( Ax)  -  (£(x  +  Ax)£(x))  - 

4  M  7T  i 

-e-\\x\!L 

(49) 

Hence,  the  mean  is  unaffected,  but  the  variance  decreases 
with  increasing  M.  Since  it  counts  the  number  of  “successful 
outcomes”  of  M  independent  trials  with  probability  ttu  C(x) 
has  a  binomial  distribution.  For  large  M  and  Trt  near  0.5,  the 
binomial  distribution  approaches  a  Gaussian  one.  Thus,  an 
advantage  of  this  description  for  global  intermittency  is  that 
it  encompasses  a  wide  range  of  behaviors,  from  two-state 
(purely  active  or  inactive)  to  Gaussian. 

The  local  spectrum  follows  from  Eq.  (16)  but  with  the 
random  C(x) ,  namely, 

C(r) 

d>G(x;/c)  =  — —  («-fl1)_ll/3_"[')/(ll/6+  v/2  ,K2a2l) 

-  y(  1 1/6  +  v/2,  K2ah)] .  (50) 

When  the  contribution  from  the  term  involving  a 2N  is  negli¬ 
gible,  /cfljS>l,  and  intrinsic  intermittency  is  absent  ( v=0), 
our  approach  reduces  to  earlier  works21’24  on  inertial- 
subrange  intermittency  where  C2  was  allowed  to  vary  with 
position.  Within  the  volume  V,  we  can  calculate  the  mean 
spectrum  <FC  by  averaging  the  right-hand  side  of  Eq.  (50) 
over  an  ensemble  of  realizations  of  C(x),  which  according  to 
Eq.  (47)  amounts  to  replacing  C(x)  with  C„tt,  . 

Figure  6  shows  two  realizations  of  QW  fields  with  glo¬ 
bal  intermittency.  The  positions  and  amplitudes  of  the  QWs 
are  actually  the  same  as  Fig.  2(b)  (a  case  with  intrinsic  inter¬ 
mittency)  but  are  then  modulated  by  a  random  C(x)  process 
with  mean  (C(x))  matching  the  original  field;  hence,  the  vari¬ 
ance  of  the  fields  is  unchanged.  The  value  of  C(x)  at  the 
center  position  b“"  is  used  to  tailor  its  amplitude  A Tan.  Fig¬ 
ure  6(a)  is  a  realization  of  C(x)l (C(x))  with  M=  1,  irl 
=  0.25,  and  L=10  m.  The  corresponding  QW  field  is  shown 
in  Fig.  6(b).  This  model  results  in  dramatic  regions  of  weak 
and  strong  activity.  Figure  6(c)  is  a  realization  of 
C(x)/(C(x))  with  M=  4,  tt^O.5,  and  L=10m.  The  corre¬ 
sponding  QW  field  is  shown  in  Fig.  6(d).  Although  the  global 
intermittency  is  much  weaker  than  Fig.  6(b),  there  are  still 
pronounced  variations  in  turbulent  activity. 


III.  COHERENCE  FUNCTION 

In  this  section,  we  derive  a  theory  for  the  coherence  of 
propagating  planar  sound  waves  based  on  the  intermittent 
QW  model  derived  in  the  previous  two  sections.  The  coher¬ 
ence  function  is  defined  here  as 
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FIG.  6.  Examples  of  synthetic  turbulence  fields  produced  by  random  QW 
ensembles  with  global  intermittency.  (a)  Random  global  intermittency  for 
C(x)/(C)  produced  by  a  Markov  process  with  77^  =0.25  and  L=10  m.  (b) 
Same  as  the  synthetic  field  in  Fig.  2(b),  except  that  the  field  is  modulated  by 
the  process  in  (a),  (c)  Random  global  intermittency  for  C(x)/(C)  produced 
by  the  average  of  four  Markov  processes  with  tt^O.5  and  L=10  m.  (d) 
Same  as  the  synthetic  field  in  Fig.  2(b),  except  that  the  field  is  modulated  by 
the  process  in  (c).  Scales  for  (a)  and  (c)  are  dimensionless;  scales  for  (b)  and 
(d)  are  in  K. 

r(x,r)  =  (p(x,r0)p*(x,r0  +  r)).  (51) 

Here,  p(x,r)  is  the  acoustic  pressure,  x  is  the  distance  from 
the  source  plane  to  the  observation  plane,  and  r0  and  r0+r 
are  points  in  the  observation  plane.  Statistical  homogeneity 
in  the  observation  plane  is  assumed,  so  that  the  coherence  is 
independent  of  r0.  As  mentioned  in  the  Introduction,  theories 
for  coherence  have  already  been  the  subject  of  much  devel¬ 
opment.  Typically,  several  important  assumptions  are  made: 
validity  of  the  parabolic  approximation,  a  long  propagation 
path  compared  to  the  size  of  the  inhomogeneities,  and 
Gaussian  statistics.  Intermittency  may  weaken  the  validity  of 
the  latter  two  approximations.  For  example,  atmospheric 
boundary-layer  thermals  and  cloud-driven  variations  in  sur¬ 
face  heating,  which  are  sources  of  intermittency,  may  extend 
over  several  kilometers.  In  this  section,  we  first  consider  con¬ 
ceptually  how  intermittency  affects  the  coherence,  and  then 
examine  simulation  results. 

A.  Local  propagation  paths 

Initially,  let  us  consider  the  coherence  when  the  propa¬ 
gation  path  length  x  is  short  compared  to  the  scale  L  of 
global  intermittency.  Paths  3-5  in  Fig.  1  illustrate  this  situa¬ 
tion.  Intrinsic  intermittency,  as  well  as  parameter  variations 
induced  by  global  intermittency,  can  impact  on  the  signal 
propagation.  For  example,  with  the  two-state  Markov  model 
for  global  intermittency  considered  in  Sec.  II  C,  the  propa¬ 
gation  path  would  be  subject  to  either  fully  active  or  inactive 
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turbulence.  This  is  important  because  we  can  consider  the 
statistics  of  the  turbulence  to  be  local  and  stable  over  the 
propagation  path. 

Previous  studies  of  intermittency  have  often  examined 
the  statistics  of  point  samples  and  line  (or  volume)  averages 
of  the  field.  The  coherence  function  actually  involves  a  non¬ 
linear  combination  of  these  types  of  statistics.  The  interpre¬ 
tation  can  be  made  clear  for  geometric  acoustics,  in  which 
the  sound  follows  ray  paths  directly  from  the  source  to  the 
receiver.  However,  interpretation  is  less  clear  when  diffrac¬ 
tion  becomes  important.  In  the  geometric-acoustics  approxi¬ 
mation,  p(x,r)=p0exp[i<j)(x,r)\,  where  p0=p  (x=0,r)  for  all 
r,  and  cf>{x)  is  a  random  phase  factor  given  by 

(f>(x,r)  =  k\  n'(x',r)dx’.  (52) 

Jo 

Here,  n'(x,r)  =  (co/c)-  1  ~-c'/cq  is  a  random  variation  in 
the  index  of  refraction.  (The  sound  speed  is  c,  c'  is  its  fluc¬ 
tuation,  and  c0  its  mean.)  For  temperature  fluctuations,  which 
in  air  are  proportional  to  the  squared  sound  speed,  n '  ~ 
—T  / 2T0  for  temperature  fluctuations.  Hence,  the  coherence 
function  is 

F(x,r)=po/exp 

(53) 


2jrk2x 

D^(x,r)  ~  — -5 —  *[1  -  J0(Kr)]<t>{K)dK.  (57) 

lo  Jo 

Here,  J0  is  the  Bessel  function  of  zero  order.  Hence,  we  have 
the  following  hierarchy:  Equation  (53)  is  valid  for  geometric 
acoustics,  the  additional  assumption  of  Gaussian  statistics 
leads  to  Eqs.  (53)  and  (55),  and  the  further  additional  as¬ 
sumption  of  a  long  propagation  path  leads  to  Eqs.  (56)  and 
(57).  One  of  the  remarkable  results  of  the  conventional 
theory  of  wave  propagation  in  random  media,  based  on 
Gaussian  statistics  and  the  Markov  approximation,  is  that  Eq. 
(54)  with  Eq.  (57)  remains  valid  in  the  parabolic  approxima¬ 
tion  even  when  geometric-acoustics  approximations  do  not. 
Put  another  way,  when  the  coherence  is  calculated  with  geo¬ 
metric  acoustics,  the  correct  general  result  is  obtained,  even 
though  the  calculation  method  does  not  correctly  capture  the 
underlying  physics  involving  diffraction  or  strong  scattering. 
Thus,  we  may  identify  limitations  to  the  validity  of  the  con¬ 
ventional  theory,  such  as  those  stemming  from  non-Gaussian 
phase  statistics  induced  by  intermittency,  by  calculating  the 
coherence  based  on  the  geometric  acoustics.  It  is  plausible 
that  the  geometric-acoustics  based  calculations  of  the  coher¬ 
ence  remain  correct  even  when  there  is  intermittency. 

Introducing  now  the  complication  that  C  may  be  a  ran¬ 
dom  function  (but  constant  along  a  particular  propagation 
path),  we  generalize  Eq.  (54)  as 

r c(x, r)  =  p\  exp[-  D^c(x,  r)/2] ,  (58) 


This  shows  that  the  coherence  function  for  geometric  acous¬ 
tics  involves  the  difference  between  two  line  averages  sepa¬ 
rated  by  a  distance  r.  A  nonlinear  function  (the  exponential) 
of  this  quantity  is  then  averaged  to  obtain  the  coherence. 
Usually,  QWs  with  size  comparable  to  r  will  most  strongly 
affect  coherence.  QWs  large  compared  to  r  are  strong  but 
tend  to  affect  both  paths  in  the  same  manner,  whereas  those 
smaller  than  r  are  relatively  weaker. 

If  the  phase  factors  r/Ax ,  r)  are  random  variables  with  a 
Gaussian  distribution,  the  formula  (exp(e))  =  exp((e2)/2) 
(where  e  is  Gaussian  with  zero  mean)  leads  to 

F (x,r)=pl  exp[-  D^x,r)l 2],  (54) 

where  Djx.r)  is  the  geometric -acoustics  phase  structure 
function,  defined  as 

D${x,r)  =  {[(f>(x, 0)  -  <fi(x,r)]2) 

k2  r  r 

=  — x  [. B(x'  -  x",0)  -  B(x'  —  x" ,r)\dx' dx" . 

2T20J0Jo 

(55) 

When  the  propagation  path  is  much  longer  than  the  size  of 
the  inhomogeneities  (x  >  a  A. 

k2x  f” 

D/.x,r)*=  — £  I  [B(x',0)  -  B(x' ,r)]dx' .  (56) 


where  and  Tc(x,r)  are  the  structure  function  and  coher¬ 
ence  associated  with  a  particular  value  of  C.  For  the  global 
intermittency  model  described  in  Sec.  II  C,  C  takes  on  dis¬ 
crete  values  Cm  following  a  binomial  distribution,  so  the  av¬ 
erage  coherence  function  is 

M 

f(x,r)  =  2  PmTcmM 

m=  1 

M 

=  Po  2  Pm  exp[-  D^cjpc,  r)/2] ,  (59) 

m=  1 

where  Pm  is  the  probability  associated  with  Cm.  If  the 
D^.c  (x,r)  are  small  (i.e.,  the  coherences  are  high  for  values 
of  C,„),  the  approximation 

exp [~D^C  {x,r)l2]~\-D(f)C  (x,r)/2  leads  to  the  follow¬ 
ing  result: 

f(x,r)  »  pi  exp[-  D/x,r)/ 2],  (60) 

where 

Drfjx,  r )  =  PJA^cJx,  r )  (61) 

is  the  average  phase  structure  function.  Since  l),b,c  “'C 
D^(C).  However,  if  the  coherence  is  low  and  the  phase 
fluctuations  are  non-Gaussian,  we  should  not  expect  Eq.  (60) 
to  hold. 

B.  Global  propagation  paths 


Using  the  spectrum,  the  preceding  result  can  alternatively  be  We  next  consider  propagation  paths  that  are  long  corn- 

written  as  pared  to  the  scale  L  of  global  intermittency,  as  illustrated  by 
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paths  1  and  2  in  Fig.  1 .  Such  paths  involve  a  varying  turbu¬ 
lence  spectrum.  A  treatment  for  this  situation  follows  from 
the  equation 


xj  /r[l- /0(/rr)]<l>(x' ; ,  (62) 


which  was  derived  in  Ref.  7  for  inhomogeneous  turbulence. 
Equation  (62)  assumes  validity  of  the  parabolic  equation  and 
the  Markov  approximation  for  the  random  medium,  and  that 
the  turbulent  fluctuations  are  Gaussian.  When  the  strength  of 
the  turbulence  varies  along  the  propagation  path  as  described 
by  Eq.  (50),  we  have  the  following  formula  for  the  coher¬ 
ence: 

f  *2 

r(x,r)  =pjj  expj  -  -^x[W(r/au\) 

[  a\*o 

C(x')cfa'|,  (63) 

where  W  is  a  deterministic  function  given  by 

_2  r°° 

W(p,s)  =  —  I  Cm-v[\-J0{fr)]y{n/6  +  v/2,es2)dZ. 
z  Jo 

(64) 


-  W(r/ai,aN/a 


l)]  f 

Jo 


exp 


D^{x,r)  r 

2x  J0 


£(x')dx' 


<exp[D0(x,r)/2], 


(68) 


The  first  inequality  follows  from  Jensen’s  inequality,  as  pre¬ 
viously  noted  in  Refs.  21  and  26.  The  second  inequality  fol¬ 
lows  from  the  argument  that  (1  lx)f£(x')dx'  cannot  be 
smaller  than  -1,  which  would  correspond  to  a  state  of  com¬ 
plete  inactivity.  Thus, 


1  T(x,  r)/pl  3*  exp[-  Djx,  r)/2] .  (69) 

Hence,  for  a  fixed  value  of  (C),  intermittency  always  in¬ 
creases  the  average  coherence. 

According  to  the  global  intermittency  model  in  Sec. 
II  C,  in  some  situations  £(x)  approaches  a  Gaussian  distribu¬ 
tion.  Actually,  since  the  integral  of  £(x)  is  effectively  the  sum 
of  many  samples  of  £(x)  when  x8>L,  it  is  even  less  restrictive 
to  consider  the  integral  as  a  Gaussian  random  variable.  As¬ 
suming  this  is  the  case,  we  obtain  the  average  coherence 
function 


r(x,r)  =pl  exp[-  £>0(x,r)/2  +  /J^(x,r)/^/8x2],  (70) 

where  I ^  is  the  following  integral: 


n*  2 

B^xj  -  x2)dxldx2  =  2 cr^Lx 

) 


1  -  -(1  -e-x/L) 
x 


(71) 


The  solution  for  this  integral  is  given  in  the  Appendix.  Note 
that  the  right-hand  side  of  Eq.  (63)  contains  a  function  (the 
integral)  of  the  random  variable  C(x).  Therefore,  to  obtain 
the  desired  expression  for  the  mean  coherence  function,  both 
sides  of  Eq.  (63)  are  averaged  over  an  ensemble  of  realiza¬ 
tions  of  this  integral.  (The  physical  meaning  of  such  averag¬ 
ing  is  discussed  in  detail  in  Ref.  21.)  By  applying  Eq.  (39), 
we  can  recast  Eq.  (63)  in  the  following  form: 


T{x,r)=pl  exp[-  D<f)(x,r)l2] 


Xexp 


Djjx.r)  r 

2x  J0 


£(x')dx' 


(65) 


where 


-  2(C)  k2x 

=  — T~T~ [W(r/ai,l)  -  !¥(>-/«!, Vai)]- 


a\To 


(66) 


The  average  coherence  function  is  thus 
f  (x,r)  =  pi  exp[-  D^(x,r)l2\ 


X  \  exp 


Dj,x,r)  r 
2x  J0 


(Xx')dx' 


(67) 


The  final  part  of  this  expression,  with  angle  brackets,  repre¬ 
sents  the  global  intermittency  effect.  It  is  bounded  as  fol¬ 
lows: 


The  final  result  for  1^  corresponds  to  the  correlation  function 
B £  given  by  Eq.  (49),  with  <r^=B {(())  =  tt{]!  M tt]  .  Equation 
(70)  provides  a  formula  for  the  coherence  function  of  a  plane 
sound  wave  propagating  in  a  turbulent  atmosphere  with  in¬ 
trinsic  and  global  intermittency.  It  can  actually  be  applied  to 
local  propagation  paths  such  that  x<SL.  Then,  e~xlL~l 
-(xl L)  +  (xIL)2I 2,  and  f^~  cj|x2;  that  is,  L  ceases  to  affect 
the  result.  The  main  limitation  of  Eq.  (70)  is  that  integrals  of 
the  function  £  along  the  propagation  path  are  assumed  to 
have  a  Gaussian  distribution. 

Figure  7  compares  coherence  calculations  from  Eq.  (70) 
with  various  combinations  of  intrinsic  and  global  intermit¬ 
tency  conditions.  The  turbulence  parameters,  based  on  the 
QW  model,  ^=0.2  m,  fli=10m,  4>l p= 0.0721,  and  ATj 
=  0.866  K.  The  path  length  is  x=  1  km,  the  frequency  680  Hz 
(wavelength  0.5  m),  and  L=500  m.  The  case  with  no  inter¬ 
mittency  (v=0  and  cr^= 0)  has  the  lowest  coherence  for  all 
separations  r.  Intrinsic  intermittency  alone  (r>=  1  /  6  and  rr; 
=  0)  increases  the  coherence  for  small  separations  (rSa{), 
whereas  global  intermittency  alone  ( v=()  and  <r!=  I  /  4)  in¬ 
creases  the  coherence  for  large  separations.  The  combined 
effect  (n=l/6  and  crt=l/4)  is  an  overall  increase  in  coher¬ 
ence;  it  is  nearly  doubled  for  large  separations. 

Figure  8  shows  the  effect  of  changing  the  length  scales 
ax  and  L.  For  each  of  the  curves,  v=\!6  and  cr^=  1/4.  One  of 
the  curves  is  for  «|  =  10  m  and  L  =  500  m,  as  in  Fig.  7.  When 
a  |  is  increased  to  20  m,  coherence  increases  for  small  sepa¬ 
rations  but  decreases  for  large  ones.  Increasing  L  to  2000  m, 
on  the  other  hand,  increases  the  coherence  for  all  separations. 
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FIG.  7.  Normalized  coherence  function  vs  the  normalized  sensor  separation 
r!ax  for  different  values  of  intrinsic  intermittency  (v)  and  global  intermit- 
tency  a ■£.  The  correlation  length  of  global  intermittency  L=1  km  and  the 
propagation  path  length  is  x=5  km. 

C.  Simulation  results 

In  this  subsection,  the  effects  of  intrinsic  and  global  in¬ 
termittency  on  the  theoretical  coherence  function  f(x,r)  are 
studied  with  simulated  random  QW  fields.  The  acoustic 
phases  are  calculated  by  numerically  integrating  Eq.  (52) 
with  the  trapezoidal  method.  The  use  of  geometric  acoustics 
greatly  simplifies  the  simulations  and  is  justifiable  for  the 
reasons  discussed  in  Sec.  II  B  2.  (Simulations  of  random 
scattering  by  QWs,  based  on  a  finite-difference  method,  are 
described  in  Ref.  37.  These  do  not  involve  the  geometric 
approximation  but  are  far  more  computationally  intensive.) 
The  simulations  are  for  essentially  the  same  situation  in  Figs. 
7  and  8,  although  with  intrinsic  intermittency  (u=l/6)  al¬ 
ways  present.  There  are  362  QW  size  classes,  and  buffer 
regions  were  used  around  the  edges  as  described  in  Sec.  II  A. 
Varying  global  intermittency  conditions  are  considered: 


FIG.  8.  Normalized  coherence  function  vs  sensor  separation  for  different 
values  of  the  outer  scale  a  {  and  the  correlation  length  of  global  intermittency 
L.  All  calculations  include  both  intrinsic  and  global  intermittency. 


FIG.  9.  Average  phase  structure  function  D^{. x,r)  at  distances  jc=12.5  m 
and  x=100  m  from  the  source,  as  a  function  of  transverse  separation  be¬ 
tween  the  observation  points.  The  outer  scale  is  10  m  and  the  correlation 
length  of  global  intermittency  is  500  m.  The  solid  line  is  a  theoretical  pre¬ 
diction.  The  other  lines  are  simulation  data:  Dashed  is  without  global  inter¬ 
mittency,  dash  dotted  is  moderate  global  intermittency,  and  dotted  line  is 
strong  global  intermittency. 

none;  M= 4,  tt^O.5,  and  L=500  m;  and  M=  1,  Tr^O^S, 
and  L=500  m.  For  the  second  of  these  cases,  cq=1/4.  We 
designate  this  case  moderate  global  intermittency.  The  third 
case,  for  which  <rl = 3 ,  is  designated  strong  global  intermit¬ 
tency.  All  statistics  were  determined  from  1024  random  real¬ 
izations. 

Figures  9  and  10  show  results  for  the  average  phase 
structure  function,  Z>^(x,r),  and  average  coherence,  f(x,r), 
for  propagation  distances  x=  12.5  m  and  x=  100  m.  At  these 
distances,  the  discussion  in  Sec.  Ill  A  regarding  local  propa¬ 
gation  paths  (x<L)  applies.  The  simulated  curves  for 
Z)^(x,r)  (Fig.  9)  are  nearly  independent  of  the  global  inter¬ 
mittency  condition,  as  they  should  be,  since  (C)  is  the  same 
for  each  case.  Theoretical  predictions  based  on  Eq.  (66) 
agree  very  well  with  the  simulations  at  x=  100  m,  but  are  too 
high  at  x=  12.5  m.  The  reason  for  this  overprediction  is  that 


FIG.  10.  (Color  online)  Normalized  coherence  function  as  a  function  of 
transverse  separation.  The  curves  correspond  to  the  same  cases  in  Fig.  9. 
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FIG.  1 1 .  Kurtosis  of  the  phase  difference  fluctuations  as  a  function  of  range 
x  from  the  source.  The  separation  between  the  observation  points  is  20  m. 
All  curves  are  from  an  analysis  of  simulation  data.  The  dashed  line  is  with¬ 
out  global  intermittency,  the  dash-dotted  line  is  moderate  global  intermit- 
tency,  and  the  dotted  line  is  strong  global  intermittency. 


FIG.  12.  Average  coherence  as  a  function  of  range  x  from  the  source.  The 
separation  between  the  observation  points  is  20  m.  The  solid  lines  are  the¬ 
oretical  prediction  without  global  intermittency  (upward  pointing  triangles) 
and  with  moderate  global  intermittency  (downward  pointing  triangles).  The 
other  lines  are  simulation  data  with  the  same  interpretations  as  in  Fig.  1 1 . 


x  —  ci\  at  this  distance,  which  means  that  the  approximation 
xi>ax  inherent  to  Eqs.  (57)  and  (62)  is  inapplicable.  As 
might  be  expected,  disagreement  between  the  simulated  av¬ 
erage  coherence  at  jv=  12.5  m  and  predictions  based  on  the 
Gaussian  approximation,  Eq.  (70)  (Fig.  10),  also  results. 
Only  one  prediction  from  Eq.  (70),  corresponding  to  no  glo¬ 
bal  intermittency,  is  actually  shown  in  Fig.  10  for  each  dis¬ 
tance.  This  is  because  the  Gaussian  predictions  are  nearly 
independent  of  the  global  intermittency  characteristics  at 
these  short  distances.  A  particularly  significant  feature  of  Fig. 
10  is  the  disagreement  between  the  prediction  and  simula¬ 
tions  at  jc=100  m  when  global  intermittency  is  present,  but 
not  when  global  intermittency  is  absent.  Since  the  average 
phase  structure  functions  are  correctly  predicted  (Fig.  9)  in 
all  cases,  the  assumption  of  Gaussian  statistics  upon  which 
Eq.  (70)  is  based  must  be  at  issue. 

The  significance  of  the  non-Gaussian  nature  of  the  glo¬ 
bal  intermittency  is  made  clearer  by  Fig.  11,  which  shows  the 
kurtosis  of  the  phase  differences,  \_4>{x ,{))-  <f>(x ,r)\  for  r 
=  20  m,  as  a  function  of  range.  The  kurtosis  for  the  case 
lacking  global  intermittency  is  nearly  3  at  all  ranges,  as  ex¬ 
pected.  For  moderate  intermittency,  the  kurtosis  decays  from 
7  near  the  source  to  about  3.5  at  x=  1  km.  For  the  strong 
intermittency,  the  kurtosis  starts  near  13  and  remains  quite 
high,  around  10,  at  x=  1  km.  One  would  expect  such  strong 
non-Gaussianity  to  affect  the  average  coherence  and,  indeed, 
as  shown  in  Fig.  12,  this  is  the  case.  The  theoretical  predic¬ 
tion  and  simulation  without  global  intermittency  are  in  good 
agreement.  For  moderate  global  intermittency,  the  simulated 
average  coherence  becomes  somewhat  higher  than  the  pre¬ 
diction  at  distances  x>  100  m.  No  prediction  is  shown  on  the 
figure  for  strong  global  intermittency,  because  the  prediction 
diverged  in  this  case.  The  simulated  coherence  is  dramati¬ 
cally  raised  by  the  strong  global  intermittency:  It  is  about 
0.55  at  1  km  versus  0.1  without  global  intermittency. 
When  there  is  such  strong  intermittency,  the  coherence  is 


dominated  by  the  many  cases  in  which  there  is  little  turbu¬ 
lence  along  the  propagation  path. 

IV.  CONCLUSION 

In  this  paper,  a  QW  model  of  temperature  fluctuations 
with  both  intrinsic  and  global  intermittency  was  developed. 
Intrinsic  intermittency  was  modeled  by  allowing  the  packing 
fraction  to  decrease  with  eddy  size.  Formulas  for  the  3D 
spectrum,  variance,  and  kurtosis  of  temperature  fluctuations 
were  subsequently  derived.  It  was  shown  that  increasing  the 
parameter  v,  which  characterizes  the  degree  of  intrinsic  in¬ 
termittency,  results  in  decreasing  variance,  but  increasing 
kurtosis,  of  temperature  fluctuations.  The  increase  in  v  also 
leads  to  steepening  of  the  inertial-subrange  spectrum  beyond 
the  familiar  -11/3  power  law. 

To  introduce  global  intermittency  into  the  QW  model, 
the  overall  turbulence  volume  was  partitioned  into  subvol¬ 
umes  in  which  intrinsic  intermittency  of  the  temperature 
fluctuations  occurs  as  described  above.  Then,  global  inter¬ 
mittency  was  modeled  by  allowing  the  temperature  ampli¬ 
tudes  AT[  of  the  largest  QWs  (eddies)  to  change  randomly 
from  one  subvolume  to  another  according  to  one  or  more 
superimposed  Markov  processes. 

The  developed  statistical  framework  for  describing  in¬ 
trinsic  and  global  intermittency  was  then  used  as  the  basis 
for  a  theory  of  the  coherence  of  a  plane  sound  wave  propa¬ 
gating  through  intermittent  temperature  fluctuations.  Calcu¬ 
lations  based  on  this  theory  show  that  intrinsic  and  global 
intermittency  both  increase  the  average  coherence.  Simula¬ 
tions  suggest  that  when  global  intermittency  is  very  strong, 
such  as  might  occur  in  stably  stratified,  night-time  condi¬ 
tions,  signal  phase  fluctuations  become  highly  non-Gaussian 
and  the  coherence  is  dominated  by  episodes  with  little  turbu¬ 
lence.  This  finding  is  important,  e.g.,  for  the  performance  of 
modern  acoustic  beam  forming  arrays,  and  should  be  exam¬ 
ined  experimentally.  It  would  also  be  valuable  to  find  ana¬ 
lytical  predictions  for  coherence  in  non-Gaussian  conditions. 
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The  QW  model  of  temperature  fluctuations  with  intrinsic 
and  global  intermittency  is  applicable  with  appropriate  pa¬ 
rameter  adjustments  to  other  intermittent  scalar  inhomogene¬ 
ities,  such  as  contaminants  in  the  atmosphere  and  density 
fluctuations  of  geological  heterogeneities  within  Earth,  and 
can  also  be  applied  to  other  types  of  wave  propagation,  such 
as  seismic  and  electromagnetic. 
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APPENDIX:  SOLUTION  FOR  THE  INTEGRAL  W 


A  solution  for  the  integral  equation  (64)  is  derived  in 
this  appendix.  We  begin  by  defining 

poo 

I(y,s,a)  =  f2a+ly(a,es2)J0(&)dl  (Al) 

J  o 

Then,  Eq.  (64)  becomes 

W(y,s)  =  (7r/2)[/(0,s,a)  -  7(y,s,a)]  (A2) 

with  a=  11/6+  v/2. 

We  now  determine  the  integral  I(y,s,a).  First,  we  use 
Eq.  (6.5.12)  from  Ref.  48  to  write  the  incomplete  gamma 
function  as  a  confluent  hypergeometric  function.  This  gives 

I(y,s,a)  =  a~ls2a  f  ^F^a'J+a;-  fs2)J0(^y)d^. 

Jo 


This  integral  can  be  solved  with  Eq.  (7.663.6)  from  Ref.  49, 
which  reads  (with  v=0) 
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j  0 


X\F\(a\b\-  \x2)J0(xy)dx 


21_T  (b) 


aay“  “e  ’  ''"'Wa/2-b+i,a/2-m 


r 
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(A3) 


Here,  is  the  Whittaker  function.  Setting  b=  1  +a  and  X 
=  s2,  we  find 


I(y,s,a)  =  2 1  -asaya-2e~r,Ss~  W_ 


■a/2,a/2-V2 


i?i-  (A4) 


We  can  use  Eq.  (13.1.33)  in  Ref.  48  to  rewrite  this  as 


I(y,s,a) 


Ur 


2  \  4 


-yW 


(A5) 


where  U  is  Rummer’s  confluent  hypergeometric  function. 
Finally,  we  use  Eq.  (13.6.28)  in  Ref.  48  to  recast  this  result 
with  an  incomplete  gamma  function  as  follows: 


Uv 


2\a-l 


2  V  4 


I{y,s,a)  =  -I  —  ]  Tll-a, — 5 
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A  practical  problem  with  this  result  for  I(y  ,s, a)  is  that  rou¬ 
tines  for  incomplete  gamma  functions  in  many  numerical 
libraries  do  not  allow  negative  arguments.  In  our  case,  with¬ 
out  intermittency  (n=0),  we  have  l-a=-5/6.  However, 
\'{a’  ,x)  [unlike  y(a’  ,x)]  is  still  a  convergent  function  even 
when  a'  <0.  This  problem  can  be  avoided  by  using  the  re¬ 
cursion  relationship 

Y{a' +  \,x)  =  a'Y{a' ,x)  +  x°  e~x .  (A7) 


Applying  this  recursion  formula  twice,  we  find 


I(y,s,a)  ■■ 


,2(a-l) 


2(a  -  1) 

1 


a -  2  4 sz 


£V-lrU  A 

a  -  2\4s2)  V  U’4s2 


(A8) 


The  argument  to  the  incomplete  gamma  function  is  now 
positive.  We  also  see  that  in  the  limit  y/s— >0, 
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Finally,  we  have 
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For  small  arguments,  T(a,z)  =  r(a)- y(a,z)  —  T(a) 
-zala.  Hence,  for  small  y/s, 
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(All) 


For  a=  1 1/6,  the  term  proportional  to  T(3 -a)  is  the  leading 
order  term,  and  we  have 
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